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On  Modeling  the  Performance  and  Reliability  of 
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Abstract 
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completion  time.  A  useful  dual  relationship,  between  the  completion  time  of  a  given 
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1.  Introduction 


We  consider  a  model  for  the  combined  evaluation  of  performance  and  reliability  of 
a  multi-mode  computer  system.  Performance  (e.g.,  throughput,  response  time, 
instruction  execution  rate)  changes  from  mode  to  mode  and  a  mode  change  occurs  in 
response  to  an  event  such  as  a  failure  or  a  repair.  The  stochastic  process  representing 
the  modes  (structure-states)  and  mode  changes  can  be  thought  of  as  a  reward  process 
by  associating  a  reward  (performance  index)  with  each  mode  [4,10].  We  can  then  study 
the  the  distribution  of  the  accumulated  reward  until  time  t  by  time  domain 
methods[l0]  or  by  transform  techniques[4,14]. 

The  authors  who  have  taken  such  a  system-oriented  view  do  not  consider  the 
effect  of  a  fault  occurring  during  the  execution  of  a  program.  A  task(job  or  program)- 
oriented  view  of  such  a  sytsem  recognizes  the  fact  that  it  is  possible  for  a  system 
failure  to  occur  before  the  completion  of  a  task  [7]  and  that  even  if  the  task  is  com¬ 
pleted.  its  completion  time  is  likely  to  be  different  from  its  execution  time  in  a  given 
mode  [3,5,12].  The  job  in  service  is  interrupted  with  each  mode  change  and  the  type  of 
service  -interruption  interaction  depends  upon  the  mode  just  entered.  For  example, 
the  occurrence  of  a  fault  during  the  execution  of  a  job  preempts  the  job  and  a  later 
system  recovery  may  allow  the  job  to  resume  from  the  point  of  interruption  (the 
preemptive-resume  (prs)  discipline)  or  the  job  may  have  to  be  repeated  from  the 
beginning.  In  the  latter  case,  the  repeated  job  may  have  the  identical  work  require¬ 
ment  as  the  original  preempted  job  (the  preemptive-repeat-identical  (pri)  discipline) 
or  a  different  work  requirement  sampled  from  the  same  distribution  (the  preemptive- 
repeat-different  (prd )  discipline). 

The  purpose  of  this  paper  is  to  develop  a  model  that  unifies  and  extends  the 
efforts  of  these  two  groups  of  researchers.  In  particular,  we  show  that  if  all  interrup¬ 
tions  are  of  the  preemptive-resume  type  then  the  completion  time  of  a  given  task  and 
the  accumulated  reward  until  a  given  time  are  dual  measures,  so  that  the  distribution 


of  one  of  them  allows  us  to  compute  the  distribution  of  the  other.  In  fact,  our  model  is 

»  i  i  \  ' 

even  more;  general  -  in  that  bottuacyilic  (closed  or  non-repairable)  and  cyclic  (open  or 

•  '  ;  ! 

repairable)  systems  are  modeled. 

Our  model  provides  an  exact  analysis  of  the  completion  time  distribution  of  a  pro- 
gram  Cob)  executing  in  a  multi-mode  system.  It  is  also  possible  to  incorporate  the 
effect  of  queueing  in  our  model.  If  the  time  spent  in  each  structure-state  is  large  com¬ 
pared  with  the  interarrival  and  processing  times  of  jobs,  then  we  can  use  steady-  state 
performance  measures  as  rewand  rates  for  each  structure-state.  Such  approximate 
decomposition  methods  have  been  considered  by  several  authors  [4,7,10,15].  If  the 
assumption  of  a  wide  separation  between  the  structure-state  holding  times  and  job  pro¬ 
cessing  times  does  not  hold,  then  a  more  complex  analysis  is  required  [1,5, 12]. 

We  develop  the  basic  modeliin  the  next  section.  In  sections  3,  4  and  5  ,  we  con¬ 
sider  the  individual  cases  where  all  structure-states  are  of  the  same  type,  that  is, 
preemptive-resume,  preemptive-repeat-identical,  or  preemptive-repeat-different, 
respectively. 

2.  77m  Basic  Modal 

Consider  a  single  server  (e.g.,  a  computer)  serving  a  single  job  (e.g.,  a  program). 
The  job  is  characterized  by  its  work  requirement,  B.  For  example,  the  work  require¬ 
ment  of  a  computer  program  can  be  measured  in  terms  of  the  number  of  instructions 
to  be  executed.  We  assume  that  B  is  a  random  variable  with  cumulative  distribution 
function  G(x)  =  P{B<x)  and  LST  \CT(s)=E(e~*B).  To  avoid  trivialities  we  assume 
G(0+)  =  0. 

The  rate  at  which  the  server  performs  work  is  assumed  to  change  with  time 

i 

according  to  the  following  model:  At  any  time  the  server  is  in  one  (and  only  one)  of  the 
n+1  states(modes)  numbered  0,t,2,.„,n.  'in  state  i  the  server  performs  work  at  rate 
r4»0,  hstscn,  work  units  per  unitltime  (e.g.,  the  instruction  execution  rate).  The  state 


3 


0  is  an  absorbing  "failure"  state,  i.e.,  once  the  server  is  in  state  0,  it  stays  there  forever 
and  the  work  rate  in  this  state  is  zero  (  r0=0).  We  allow  absorbing  non-failure  states 

among  the  states  1 . n  with  reward  rates  greater  than  zero  so  that  if  the  server  enters 

such  a  state,  the  job  will  eventually  complete.  Let  Z{t)  be  the  state  of  the  server  at 
time  t.  \Z{t),t  ^  0}  is  called  the  structure-state  process.  We  shall  assume  that  the 
structure-state  process  is  a  stochastic  process  with  piecewise  constant  paths  with  fin¬ 
ite  number  of  jumps  in  finite  intervals  of  time.  Furthermore,  the  structure-state  pro¬ 
cess  is  assumed  to  be  independent  of  the  work  requirement  B  of  the  job. 

The  states  i  =  1,2,. .,n  are  classified  as  (i)  prs :  preemptive-resume,  (ii)  pri: 
preemptive-repeat-identical  or  (iii)  prd:  preemptive-repeat-different. 

The  following  quantities  have  been  analyzed  before  in  the  literature  for  some  spe¬ 
cial  \Z(t),t  >  Oj  processes: 

I.  The  job  completion  time  (T(x)):  defined  to  be  the  total  time  the  server  takes  to  com¬ 
plete  a  job  that  requires  x  units  of  work.  T  denotes  the  unconditional  completion  time 
of  a  job  that  requires  a  random  amount  of  work,  say  B  units.  Gaver[5]  studied  the  dis¬ 
tribution  of  the  r.v.  T  for  a  system  subject  to  one  type  of  failure  and  repair,  in  which 
the  operating  state  is  Markovian  and  the  failure  state  is  semi-Markovian.  Nicola[l2] 
extended  Gaver’s  model  to  allow  for  mixed  types  of  failures  and  repairs.  Castillo  and 
Siewiorek  [3]  condsidered  a  system  with  two  types  of  failures  in  which  the  preemptive- 
repeat  type  failure  could  occur  during  the  repair-time  of  the  preemptive-resume  type 
failure. 

II.  The  probability  of  dynamic  failure  ( rj ):  defined  to  be  the  probability  that  the  sys¬ 
tem  fails  before  the  job  is  completed,  Le.  the  server  enters  state  0  before  completing 
B  units  of  work[7]. 

III.  77ie  cumulative  reward  upto  time  t  (Y(t)):  defined  to  be  the  total  amount  of  work 
done  by  the  system  up  to  time  t .  Y  is  the  total  accumulated  work  during  the  system’s 
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lifetime;  it  is  the  limit  of  Y(t )  as  t  -»  <*>.  The  r.v.  Y(t )  was  first  studied  by  Puri  [  14]  for 
Markovian  Z(t)  processes.  Meyer[lO]  and  Donatiello  and  lyer[4]  studied  the  distribu¬ 
tion  of  y(f )  for  an  acyclic  Markovian  Z(t)  process.  Beaudry  [2]  studied  the  r.v.  Y  for  a 
Markovian  Z(t)  process,  while  Osaki  and  Nishio  [13]  studied  the  r.v.  Y  for  a  semi- 
Markovian  Z(t)  process. 

To  present  a  unifying  view  of  the  quantities  defined  above,  we  introduce  the  cumu¬ 
lative  measure,  W(t),  defined  as  follows:  Suppose  that  at  time  t  =  0  the  server  starts 
processing  a  job  with  infinite  work  requirement.  Jf(f)  is  the  amount  of  useful  work 
completed  by  the  server  until  time  t  (thus,  excluding  the  work  done  prior  to  the  last 
visit  to  a  pri  or  apnf  state).  The  following  properties  of  the  cumulative  measure,  W(t), 
are  immediately  obvious: 

(i)  *(0)  =  0. 

(U)  Z(t)  -i  =>  dW(t)/dt  =  rt, 

(iii)  If  there  is  a  transition  in  the  structure-state  process  at  time  t  and  Z(t  +)  =  i, 

then  W(t+)  =  0  if  i  is  a  pri  or  a  prd  state  and  tf(f+)  =  tf(f-)  if  i  is  a  prs  state. 

Typical  sample  paths  of  the  structure-state  process  and  the  cumulative  measure, 
are  shown  in  figure  1,  for  the  following  case:  Set  of  states  =  ]0, 1,2,3],  state  1  is 
prs  with  =  1,  states  2  and  3  are  pri  or  prd  with  r8  =  2  and  r8  =  0,  state  0  is  the 
absorbing  failure  state. 

The  following  theorem  shows  how  the  quantities  T,  r),  Y(t)  and  Y  can  be  related  to 
each  other  via  the  cumulative  measure  ,  W(t ). 

Theorem  1. 

(i)  rsminjf  Jfc  0:  Jf(f)  =  flj. 

(ii)  The  dynamic  failure  probability,  fj  ~  P(T  =  <*>). 

(iii)  If  all  states  are  prs,  then 


■r'7-sT.. 


-  •  •  *.%’•  j*  * 


v-./.  • .  -v 


P(Y(t)<x)  =  l-P(T(x)<t) 
and 


P(Y*x)  =  l-P(T(z)  <«). 

Proof:  (i)  Let  T  be  the  job  completion  time.  It  is  clear  that 
\T  >t\ *±\W(u)<Bt  for  all  OsJusJf  J, 

since  W(u)  represents  the  useful  work  done  upto  time  ti.  As  W(t)  has  piecewise  con 
tinuous  paths  with  only  downward  jumps,  T  is  given  by  (i). 

(ii)  It  is  clear  that 

(  Dynamic  Failure  J  «=»  \  system  fails  before  job  completion  J 
<*»\W(t)<B  for  all  t  is  Oj 

Hence  rj  =  P(T  =  «). 

(iii)  Let  T(x)  -  minjf  >  0:W(t)  -  x\.  If  all  states  areprs,  then 

{ Y{t )  >  x  l  «=►  ( fT  ( t )  >  x  l  \  T(x )  <  t  i . 

Hence 

P(Y(t)  >  x)~  P(T(x)  <  t).  Q.E.D. 

It  is  apparent  from  the  above  theorem  that 

F^minlt  =  (2.1 

is  the  unifying  random  variable.  This  paper  is  devoted  to  the  study  of  this  random  vari 
able.  Define  the  following  distribution  functions: 

Fi(t,x)  -  P(T  <  t  \B  =  x,  Z(0)  =  i),  *2*0,  lsisn, 

F(t,x)  =  P(T<t\B  =  x),  x>0, 

Fx(t)  -  PIT  <  t  \Z(0)  -  i),  Kisin, 


.•  sxs.  »  -v  * 


‘  \r\  ' 


F(t)  =  P(T*zt) 

and  the  corresponding  LSTs  (Laplace  Stieitjes  Transforms), 

/l~(s jx)  =  ^(e-*r|  B  =  x,Z(  0)  =  i),  iiO,Uisn,  (2.2) 

F“(s.z)  =E(e-*T\B  =x),  xfcO,  (2.3) 

/T(s)  =  ^(e  ~*T  |  Z(  0)  =  i).  lsiin,  (2.4) 

F~(s)  =  E(e~*T).  (2.5) 

From  the  independence  of  fi?(f ).*  &  Oj  and  5  it  follows  that 

r(s,x)  =  £  Fr(s,z)P(Z(0)  =i),  z  >  0.  (2.6) 

4*1 

/T(s)  =  J  Fi~(s  jc)dG(x),  1  Si<n,  (2.7) 

0 

/~(s)  =  £  /VW(Z(0)  =  i).  (2.B) 

4*1 

From  equations  (2.6)  -  (2.B)  it  is  clear  that  the  conditional  LSTs  F?(s,x)  are  of  central 
importance  to  the  analysis  of  T.  In  order  to  obtain  explicit  formulae  for  FC(s  ,x )  it  is 
necessary  to  make  some  further  assumptions  about  the  structure-state  process.  In 
the  remaining  paper  we  make  the  assumption  that  \Z(t  ),t  >.  Oj  is  a  time  homogeneous 
continuous  time  Markov  chain  (CTMC).  The  results  derived  here  can  be  extended  in  a 
straight  forward  manner  to  the  case  when  the  structure-state  process  is  assumed  to  be 
semi-  Markov.  Let  gy,  1  <i  *  )  in,  be  infinitesimal  transition  rate  from  state  i  to  j 
and  qio  be  the  absorbing  failure  rate  from  state  i.  Let  Q  =  [gy],  1  s  i,  j  «n,  be  the  n 

by  n  generator  matrix  where  g<  =  £  gy  =  -g« .  Note  that  row  sums  of  Q  are  <■  0.  We 

f=o 

/»•< 

mention  one  property  of  the  CTMC  for  future  reference.  Define 

H  =  min(f  Js  0;Z(f )  *  Z(0)|  (2.9) 

as  the  holding  (or  sojourn)  time  in  the  initial  state.  Then  we  have 

(  ~)  denotes  LST, the  Laplace  transform  of  a  probability  density  function. 


"v  "  •*  *.  a  ".a  *  -a  0  m  0  h  "  >■  *  »  •  — m  '  A  *  • V  m  *  *  •_>  k  ■  «  "  ■  *»***'•  *’*»*_•  *  •  "  » 
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P(HHx,Z(H+)=j\Z(0)=i)  =  ^-(l-e-*').  (i  *  j).  (2.10) 

9t 

In  the  next  section  we  treat  the  case  where  all  states  i  =  1,2,..ti  are  preemptive- 
resume  (prs)  and  in  sections  4-  and  5  we  consider  the  case  where  all  states  are 
preemptive-repeat  (pri  and  prd,  respectively).  The  mixed  cases  where  some  states  are 
prs  and  some  are  pri  or  prd  have  been  studied  in  [8]. 


3.  The  Preemptive-resume  Case 

In  this  section  we  assume  that  the  states  1,2, ...n  are  all  preemptive-resume 
states.  Note  that  state  0  does  not  have  to  be  classified  since  it  is  a  failure  state. 
Theorem  2  below  gives  a  method  of  computing  the  conditional  LST s  defined  by  equa¬ 
tion  (2.2).  First,  some  notation: 


FC\s.u)  =fe~**  Fi~(s.x)dx  ,  1  £  i  <  n. 

0 

r-(s.ti)  =  [/T(s.«).  F;\s.u) . F- *(s,u)]r, 

R  =  diag[rl%rz . r„]  . 

X  =  iTi.rg . r„]r, 


where  the  superscript  T denotes  transpose. 


Theorem  8.  FC\s,u)  .for  1  <  i  •£■  n  ,  is  given  by 

Fi\s,u)  -  — — +  £  ■  --V— — FfXs.u),  lsiin, 

*  v  s  +r4u  4-g*  /^s+nu+gt  • 

in 

Proof:  Conditioning  on  the  sojor~n  time  H  in  the  initial  state  we  get 

if/iisx/rt 


E(e~*T\H=h,B-x,Z(0)=i)  = 


e^f]  ^-Ff(s,x-rih),  if  h  <x/rt 
in 


(3.1) 

(3.2) 

(3.3) 
(3-4) 


(3.5) 


(•)  denote*  the  Laplace  transform  of  a  function 


Unconditioning  yields 


FTisjc)  =  /  E(e~tT\H  =  h,B  =  x,Z( 0)  =  dA 

0 

=  „>*«>'»  fr(s.x-r,h)dh 

}S  0 

Multiplying  both  sides  by  e  and  integrating  we  get  equation  (3.5).  Q.E.D. 

Equation  (3.5)  can  be  put  in  a  matrix  form  as  follows: 

[si +  uR  -Q)r'(s.u)=Z. 

where  I  is  the  identity  matrix.  As  it  is  well  known  that  [si  +  uR  -  9]  is  invertible,  we 
get 

£~*(s,u)=[s/  +  u/?  -  C]_,H  •  (3.5a) 

A  direct  inversion  with  respect  to  s  yields 
df  £*(f  .it)  =  e^Q~uR^t  r.  . 

After  integration  and  some  manipulations,  we  get 

£‘(f,u)=  l-[/- .<•-««*] j,  .  (3.6) 

l JL 

We  now  describe  how  we  can  use  the  above  theorem  to  compute  Ff(s,x).  Using 
Cramer's  rule  we  can  write 

Fi\s,u)  -  At (s.u)/  C{s,u) 

where  C(s,u)  =  det[sl  +  uR  —  9]  and  Ai(s,u)  are  appropriate  n  by  n  subdeterminants 
of  the  augmented  matrix  [s/  +  uR  -  ^;xl-  It  is  obvious  that  both  Ai(s,u)  and  C(s,u) 
are  polynomials  in  s  and  u.  Hence  one  can  use  partial  fractions  to  invert  Fi"\s,u)  with 
respect  to  u.  Let  d  =  |(i:rt  >  0j  | .  i.e.  d  is  the  number  of  states  in  which  work  rate  is 
positive.  Then  C(s,u)  is  a  d-degree  polynomial  in  it  for  a  fixed  value  of  s.  Let 
—Ui(s) . — Urf(s)  be  the  roots  of  C(s,u).  In  the  special  case  when  these  roots  are  dis- 


tinct,  we  can  write 


FCXs.u)  =  £ 
1=1 


Ai(s) 

U+Uj(s)' 


l^i^n  , 


where 

s) = w^r(u + Ui(s))'  1  * d 

Inverting  with  respect  to  u.  we  get 

FC(sj:)  =  £  v4y(s)  .  Ui^n. 

l=i 

Hence  from  equation  (2.4) 


(3.7) 


(3.8) 


(3.9) 


K(s)  =  £  i4y(s)  CT(uj(s))  .  lsisn,  (3.10) 

1=i 

(recall  that  CT(s)  =  f  e~**dG(x)),  and 
o 

r-(s)  =  £  [£  rr4>ly(s)]  G~(uj(s)).  (3.11) 

1=i  i=i 

where  rq  =  P(Z( 0)  =  i),  1  5i«n. 

It  is  interesting  to  note  that  the  LST  of  T  for  a  given  s  is  simply  a  linear  combina¬ 
tion  of  the  LST  of  B  evaluated  at  itj(s) . lirf(s). 

Now,  assuming  that  state  0  is  reachable  from  every  other  state,  the  probability  of 
dynamic  failure  can  be  computed  easily  from  Theorem  1  as 


T)  =  p(T  =  «)  =  l  _  lim  F~(s).  (3.12) 

« -*o 

The  following  corollary  indicates  how  the  LST  of  the  cumulative  reward  Y(t ),  for  a 
given  t,  can  be  obtained  from  the  F<~*(s,u)  functions. 


OrroUary  1.  For  a  given  t  St  0,  let  Y(t)  be  the  cumulative  reward  upto  time  t .  Let 


Yi(x.t)  =  P(Y(t)^x  I  Z(0)  =  i), 


10 

Y?(u,t)  =  I  Z( 0)  =  i) 

and 

Y?*(u,s)=f  e~«  Y?(u,t)dt  . 
o 

Then 

^•(u.s)  =  ^-(1  -uFC*(s,u))  ,  lsisn,  (3.13) 

Proof:  Part  (iii)  of  Hieorem  1  implies  that 

P(Y(t)  <  x  |  Z( 0)  =  i)  =  P(T(x)  >  f  |  Z{ 0)  =  i). 

Now, 


YT*(u,s)  =f  z(0)=i)dt 

0 

=  /e-7  8““d,^(Ir(O«*l^(0)=i)df 
/~!0  »*0 

=  /  e-«7  «",‘4^(no=s*i^(o)=i)df  • 

r=0  r=0 

=  7  e-«dx[7  e-‘[l-P(r(x)^f|Z(0)  =  i)]ctt] 

«=0  1=0 

=  -  7  *““<4  FC(s,x)/s  =  [1  -u  /V*(s,u)]/s.  Q.E.D. 

0*0 

Using  equation  (3.5a),  we  can  write  in  a  matrix  form 

iT*(u,s)  =  [s/ +  u/?  -  9]"1  J. 
with 

JT*(«.s)  =  [lT*(«.s).17*(u.s) . !£*(u.s)]r. 

A  direct  inversion  with  respect  to  s  yields 
JT(u,t)  =eW -**)*&  . 


(3.13a) 


(3.14) 


SS 
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We  end  this  section  with  a  simple  example. 

Example  3. 1.  The  switching  server 

Consider  a  system  that  operates  in  two  modes  each  with  a  different  work  rate,  say 
T\  and  re  for  modes  "1"  and  "Z\  respectively.  The  system  switches  between  the  two 
modes  according  to  a  Poisson  process  at  different  rates,  say  X  and  pi  from  modes  "1" 
and  "2",  respectively.  A  total  system  failure  may  occur  at  any  mode  of  operation  at  dif¬ 
ferent  rates,  say  Xo  and  pio  for  modes  "1“  and  "2”,  respectively.  The  CTMC  representing 
the  switching  server  is  shown  in  figure  2.  In  the  case  where  a  total  system  failure  may 
not  occur,  i.e.  X<j  =  piQ  =  0  and  if  r8  =  0  then  the  switching  server  model  reduces  to  the 
completion  time  model  of  job  execution  in  a  system  subject  to  breakdowns  and  repairs 
considered  by  Gaver  [5]. 


In  this  example  we  consider  the  case  in  which  both  states  1  and  2  are  of  the 
preemptive-resume  type.  We  note  that  if  we  set  pi  -  0  in  this  example  we  obtain  the 
reward  model  of  a  two  processor  system  considered  by  Meyer  [10],  In  our  example  the 
Q  matrix  is 


I 


* 


f7(s,x)  =  A,j(s)«xp(— u,(s)*)  +  A18(s)«xp(-u8(s)*), 


\z 


/,£,(s.x>=i4tl(s)837>(-u1(s)jc)  +  Azi(s)exp(~U2(s)x) 


where 

1 


fttj(s)  =  tr1(s+/i')+r2(s  +X,)+V(r1(s  +n’)-r^s  +X’))z+4X^r  ,r2]/  (*r  ,r  j) 

i  1  I 

=  [r,(s  +/*’)+r2(s  +X’)-V(r  ,(s  +p,')-r2(s  +X’))z+4\ur  ,r2]/  (Sr  ,rjJ) 

Au(*)  +  -«i(s))»‘i*,2] 

i4l8(s)  =  (r,(jf +p’)  +  r2X  -  T,r*u2(s))/[(u,(s)  -  u2(s))r,r2] 

Auk)  =  lre(s+\')  +  -rjrgu,(s))/[(u2(s)  -  u^sJJrjTa] 

XasOO  =  (r2(s+X')  +  r ,/i  -  r  jrgu^s  ))/  [(u  ,(s  )  -  ix2(s))r,r2] 

Then 

F"(s)  =  [Mjj^J+Trs^zj^^^u^sJJ+frr^^sJ+rrz^sJjGr^s ))  .. 
And  17,  the  probability  of  dynamic  failure,  is  given  by  1  -  F~( 0). 

From  corollary  1,  we  have 

>?•(«.*&  =  j-[i-«Fr#(s,«)] 

_  (a  +X’)(s  +fX)  +  reu(s+Xc)  ~  V* 

«[(*+X'+rjii)(s  +/4*+r2u)-^i] 

=  Pio(u)  f  ^n(u)  |  BjZ(u) 

S  S+Sj(lf)  s+s2(u) 


where 


*i(u)  ~  +  riu  +  M'  +  rzw)  +  V(X'+r1ii-/i'-r2u)*  +  4X/i)] 


*eM  =  g^X'  +  r,,u  +  +  t2u)  -  y/(\’+rlu-/u,'-rt u)z  +  4\^)] 


Ka(u)se(u) 


B„(u)  = 


_  K*'"«lM)(M'-»t(M))  +  r2u(X0-s1(u))  -\m 


c,(u)[s,(u)-sB(u) 


,  .  KA'-f^u^Qi'-SgCu))  +  r2u(Xo-s2(u))  -  X/x 

^  ’  ~  *£<«)[S8(“)-Sl(w)] 


Inverting  with  respect  to  s ,  yields 

YUu.t)  =  Bt  0(u)  +  Bll(u)e~,'{u)t  +  B^u)*'**'1* . 

In  a  similar  manner  we  can  compute  Yg(u,t).  We  note  that  the  above  LSTs  can  be 
inverted  in  this  case  to  obtain  the  distribution  function  of  Y(t)  as  an  infinite  sum  of 
Bessel  functions  owing  to  the  occurrences  of  radicals  in  the  expresions  of  Si(u)  and 
se(u).  However,  in  the  case  that  \i- 0  (as  considered  by  Meyer),  the  radicals  disappear 
and  the  inversion  is  relatively  easy  (as  derived  in  [4]  for  arbitrary  number  of  proces¬ 
sors). 


4.  77ie  Pre emptive-repe at^identic al  Case. 

In  this  section  we  assume  that  all  states  are  preemptive  repeat-identical.  The 
m«un  result  is  given  in  the  following: 


Theorem  3.  The  conditional  LSTs  Ff(s,x),  1 «  i  <  n  as  defined  in  equation  (2.2)  satisfy 
the  following  simultaneous  equations: 


*-(«.*)  =  +  S  i  *  t  *  (4J) 


Proof  ■  Conditioning  on  the  holding  time  H  in  the  initial  state  we  have 


E{a-'T\H=h.B=x,Z{0)=i)  = 


t  ,r<  .ifhfcx/r* 

^-Fr(s,x),  if  h  <*/rt 
>- 1  7< 
in 


Unoonditioning  yields  equation  (4.1).  Q.E.D. 

Solving  equations  (4.1)  we  get  /V'(s.x),  for  ltfi  £n.  Then  equations  (2.7)  and 
(2.8)  can  be  used  to  compute  F~(s).  Finally,  rj  =  1  -  F“(0). 


Example  4. 1.  Consider  the  switching  server  of  example  3. 1,  except  now  we  assume  that 


states  1  and  2  are  preemptive-repeat-identical. 

Equations  (4.1)  become: 

(s+\>)F:(s,z)  =  +  X(1  -e-l,*X)*/r')FZ(s,z) 

(s+fi‘)Fz(s,z)  =  +  Md-e-{,^/r«)Fr(s.x). 

Solving  the  above  equations  we  get 

/7(s,x)  =  (s+/Li')[s+X')a  +  X(l-a)6]/A 
Fz  (s,x)  =  (s +X)[s +/x')b  +  /u(l-6  )a]/A 

where  a  =  eajj^s+X'Jx/f!),  b  =  ezp(-~(s+p')z/r  2)  and 

A  =  (s+X')(s+/*')  —  X/i(l-a)(l-6)  .  Fi~(s),  for  i  =  1,2  and  F~(s)  can  be  obtained  from 
equations  (2.7)  and  (2.8). 

5.  The  Preemptive-repeat-diff event  Case 

Here  we  consider  the  case,  where  all  structure-states  of  the  process  are 
preemptive-repeat-different  (prd). 

The  following  theorem  holds 

Thmovem  4.  The  LST s  F?(s),  for  1  <  i  <  n,  as  defined  in  equation  (2.4)  satisfy  the  fol¬ 
lowing  simultaneous  equations 

/TOO  =  C~((s+9<)/i\)  +  J  ^j^-y{l -tr((s+g<)/r<)]F/'(s),  1  *i*n.  (51) 

Note  that  when  rt  -*  0,  G~((s+g<)/r<)  -»  0,  since  G  (0+)  =  0  and  hence  LimG~(s)  -*  0. 
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Proof:  Conditioning  on  the  work  requirement  B  of  the  job  to  be  executed  and  on  the 
holding  time  H  in  the  Initial  state  we  get 


E(M-*T\B=z.H=h.Z(0)=i)= 


if  A  fex/r* 
if  h  <  x/rt 


j-i  7* 


Note  that  if  a  structure  state  transition  occurs  before  the  job  is  completed  then  a  dif¬ 
ferent  job  with  independent  and  identical  distribution  is  restarted. 

Now,  unconditioning  on  B  (the  job’s  work  requirement)  yields 


E(fT\H=h,Z(0)=i)  =  f  #“**/ridG(x)  +  ?  (s)dG(x) 

-  fm  1  9i 

J+i 


«*0 

Unconditioning  on  H  (the  holding  time  in  the  initial  state),  yields  equation  (5.1). 
Q.E.D. 


Solving  equations  (5.1)  we  get  Fi~(s),  for  lcisn.  Equation  (2.8)  can  be  used  to  get 
F"(s).  The  dynamic  failure  probability  (rj)  follows  immediately 

rj  =  p(7’=«)=l-r'(0). 

Note  that  the  preemptive-repeat-different  case  with  a  constant  (or  deterministic) 
work  requirement  of  a  job  (B=x)  corresponds  to  the  preemptive-repeat-identical  case. 

Example  5. 1.  Again  we  consider  the  switching  server  of  example  3. 1  with  the  states  1 
and  2  being  preemptive-  repeat-different.  From  equations  (5.1)  we  have 


F,~(s)  =  CT((s+  A’)/r,)  +  ^y{l -(T((i +*)/*•,)]*?•(») 

rg~(s)  =  C~((s+M')/re)  +  1  -G~((s  +M')/  /  r  e)]Ff  (s ) 

It  follows  that 

CT((n  +A')/rl)+(-~— )(l~C~((s  +A')/r,))G~((s  +/i’)/re) 

*7(»)  =  - 7 - - 

[Mj ^(^(1“G~((s+X’)/rI))(1-fr((s+M’)/re))] 
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f;(s)  = 


cr«s  +fi‘)/rt)H-£zHl-Cr((S  ^)/rg))(T((s  +\')/rl) 

[H^Q(1'<r((s+v)/rl))(1"(r((!+'1’)/r8))1 


F~(s)  can  be  obtained  from  equation  (2.8). 


6.  Conclusions  and  Extensions 

We  have  developed  a  unified  model  for  the  combined  evaluation  of  performance 
and  reliability  of  multi-mode  computer  systems.  This  allows  us  to  compute  both 
system-oriented  measures  (such  as  the  accumulated  reward)  and  task-oriented  meas¬ 
ures  (such  as  the  completion  time  and  the  dynamic  failure  probability)  from  a  single 
model.  We  model  preemptive-resume  and  preemptive-repeat  interactions  between 
task  execution  and  mode  change  (failure/repair)  events.  It  is  clearly  of  interest  to 
allow  mixed  preemptive-resume  and  preemptive-repeat  interactions  in  the  same 
model.  Tins  and  other  extensions  have  been  studied  and  reported  recently  [8].  The 
techniques  developed  in  this  paper  can  be  exetended  to  the  case  where  the  structure- 
state  process  is  a  semi-Markov  process. 
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Hence,  using  eq.  (3.9)  we  get 

F?(s,z)  =  ^n(s)eip(-Ui(s)i)  +  X^sJeapC-UgCsJz), 
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Solving  equations  (4.1)  we  get  Fi(s,x),  for  lsiift.  Then  equations  (2.7)  and 
(2.0)  can  be  used  to  compute  F~(s  ).  Finally,  rj  =  1  -  F~(0). 

Example  4.1.  Consider  the  switching  server  of  example  3.1,  except  now  we  assume  that 


* 


V  O  VAVs'aV.V/  v  O  O  V'  c  Oa.»  aJ 


